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Abstract 

Most of the thermodynamic analysis of ovens are focused on efficiency, 
but they need to behaves under real-hfe conditions, then the effectiveness 
of the ovens plays a crucial role in their design. In this paper we present 
a thermodynamical model able to describe the temperature evolution in 
ovens, furnaces or kilns to harden, burn or dry different products and 
which provides a methodology to design these heating devices. We use 
the required temperature evolution for each product and process as main 
ingredient in the methodology and procedure to design ovens and we place 
in the right role the efficiency criteria. We use global energy balance 
equation for the oven under transient situation as the thermodynamic 
starting point for developing the model. Our approach is able to consider 
different configurations for these heating devices, or recirculating or open 
situations, etc. 

1 Introduction 

The brick kiln was a major advance in ancient technology because it provided 
a stronger brick than the primitive sun-dried product. Modern kilns or ovens 
are used in ceramics to fire clay and porcelain objects, in metallurgy for roast- 
ing iron cores, for burning lime and dolomite, and in making Portland cement 
between other applications. This kind of ovens are an essential part of the man- 
ufacture of all ceramics, which, by definition, require heat treatment, often at 
high temperature. During this process, chemical and physical reactions occur 
which cause the material to be permanently altered. In the case of pottery, 
clay materials are shaped, dried and then fired in a oven. Although there is 
an experience of some several thousands years building ovens or furnaces and 
we can find some high technology apparatus, very detailed numerical analysis 
(for example see jSj), two phase flow study (for instance ^) or efficiency based 
analysis [5]; up to our knowledge there is no a simple model that allows us to 
calculate the physical properties and design an oven. 




Figure 1: Schematic oven and its main parts: 1. Chamber, 2. Thermal insulator, 
3. Energy source, 4. Material to be burn. 
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Figure 2: Cooking process with three major phases in the heat treatment, in- 
creasing temperature, cooking period and coohng period. 



In figure [U we present a scheme of our conceptual oven composed by an insu- 
lated chamber with a thermal energy source in which the material to be burned 
can be introduced. The heat can be obtained by different energy sources those 
can be inside or outside the chamber, for the purpose of this paper the energy 
source is not important. The schematic oven is composed by a cooking chamber 
which is thermally insulated with a high thermal resistivity material; inside this 
chamber the material to heat treated is placed. Both material and chamber are 
heated simultaneously during oven operation with internal or external energy 
source. 

The aim of this paper is not only to describe the thermal behavior of a 
heating device, but to establish a methodology for ovens design considering the 
temperature evolution for different products and process. 

A typical temperature curve of cooking [TT] is given in Fig. [21 where we 
can see that the process starts at i = with temperature T^. If the process 
consists in a constant heat flow i?o,i as indicated in Ilia during a time tc, the 
internal temperature will follow a curve as la up to reach To then a maturity 
period follows, as shown in part lb- Of course to follow this plot the heat 
flow needs to be diminished. Paschkis and Persson |Tl] proposed a slow change 
in heat flow as indicated in Illb p. In order to obtain simplest model we will set 

this heat flow as constant, i?o,/ (-E'o.j > ^o,/) as plotted in curve Ilh- After 
the cooking process the cooling part follows with E = 0, curve ///c, then the 
temperature may be described by a curve like /c, at the end of cooking the 
thermal energy flow to the oven is canceled (curve I lie)- With this the cooling 
of the oven and the products start (curve Ic), and the temperature ends in the 
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temperature at time te- In the same figure, it is shown that if the operation 
temperature To was not specified (in other words To < T^ax), and the energy 
is introduced i?o,i to the oven as shown on the hues Ilia and Hid, the thermal 
process starts foUowing the temperature curve la until time tc, and later, it will 
continue the curve II. We can see that the curve // has as asymptotic limit 
to the line of constant temperature T^ax, which is reached only when t ^ oo. 
Of course, the graph of the figure [5] has been observed in many experimental 
works in the firing of miscellaneous products in ovens furnaces or kilns, when 
the average temperature is measured, and has obtained analytically data to 
interpolate; however, there is not a well-founded theoretical model. 

We must emphasize the efforts to shape the curve of temperature with simple 
models. In 1961, Trinks and Mawhinney [17j . modeled the temperature curve 
of the load when the oven is preheated and the load is introduced at room 
temperature. Then in 1994, Astirraga [5] presented a simple model for the 
thermal development of an electric oven, and he founded the differential equation 
that describes the process of warming, however, the initial conditions of the 
thermal process were not well established. Therefore, the found solution does 
not correspond to the baking process, and this only describe qualitatively the 
process. In 2002, Abraham and Sparrow [1], described the thermal behavior 
of a load at room temperature introduced on an electric preheated oven; the 
analytic relation found has similarity with Trinks and Mawhinney's work |17| . 
In 2005, Tapia and del Rfo [16] characterized and model theoretically a "solar 
cooker" finding that the analytic representation of the temperature evolution 
inside the cooking chamber of the cooker depends on the thermal resistance of 
the oven, its heat capacity, and the characteristic parameters of the process. 
They also showed that cooking of food had a thermal development according 
to the predictions of theoretical models. In 2007 Schwarzer and Vieira [14] 
developed a theoretical model, to evaluate the thermal development in a pot 
heated by solar energy, however, they apply the relation for two different cases, 
the first one for cooking in a pot exposed to the atmosphere, and that receives 
the solar flux concentrated by the reflect of a mirror, and the second one to a 
solar oven box. We should emphasize that, according to their approach, their 
theoretical model, can only describe the balance of energy in the pot exposed 
to the environment, and not to the pot inside the oven. 

The goals of this paper are: develop a methodology to design ovens based 
on a simple thermodynamic model that describes the temperature evolution 
according with the heating process indicated in figure. [2] 

2 Modeling typical heating curves in ovens 

In this section we present a theoretical model for describing the temperature 
evolution inside of an oven in terms of some parameters that allow to us to 
design ovens or furnaces according to specific cooking requirements. 
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2.1 Describing temperature evolution. 

We wish to describe the evolution of a representative temperature inside the 
cooking chamber of a oven. Here it is important to mention that a detailed 
description of the temperature field may be important for some specific processes 
or for big cooking chambers where stratified temperature conditions could be 
found. However in this first approach, we can assume that we are selecting 
a good place to measure the temperature and it is representative of the whole 
temperature field, with this representative temperature we can construct an one 
temperature model for describing the cooking process as we have done above. 
In order to do the one temperature model we will use a global balance of energy 
using the cooking chamber as the control volume. Due to the fact that, we try 
to reproduce the cooking curve, shown in figure [2l we will consider the next 
statements: 

1. We consider that all objects in the chamber are under thermal equilibrium 
condition. 

2. The materials that compose the oven and the cooking material have the 
constant physical properties, i.e. they do not depend on temperature. 

3. The heat capacity Cj = CpjmjUjj of insulators and supports of the oven 
are considered in the balance energy equation, where Cpj, mj are the 
specific heat and the mass of the component j in the oven. 

4. The change of temperature T with respect to t has the following require- 
ments: ^ is maximum at t = 0, it is a non increasing function of t and 

— > when t — > cx). 

at 

5. We neglect the energy of the chemical or phase transformations. 

6. We assume that no work enters or goes out from the furnace 

1. We consider that cavity and both air and body surfaces places in it are in 
thermal equilibrium. 

2. The physical properties, as thermal conductivity and specific heat of oven 
components and cooking materials, are constant. 

3. The temperature change in the inner part of the oven satisfies the following 
restrictions: 

• £ is small in the heating process. 

• e(t) is maximum at t = 0, i.e., it is a decreasing function on t. 

• e{t) when t ^ oo. 

4. We consider the following approximations for the cooling part. 
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• I e I is small aX t ~ tc ~\- tm- 

• I e(t) I is maximum at t = tc + tm, i-e., it is a decreasing function on 
t. 

• I s(t) |-> when < — s- oo. 

5. We are neglecting all the energy involved in chemical changes 

Following these ideas, the energy balance equation of the oven can be ex- 
pressed as: p^-[20]: 

Qe~Qs^^, (1) 

where Qe is the input heat, Qs is the heat diffused into the environment, Eac 
is the accumulated energy of the furnace. 

Here we can use that the accumulated energy can be expressed as 

^ = C.^ (2) 
dt dt ^ ' 

where Ct = Cp,jmjUij is the effective heat capacity of the oven. Here we 
are considering the supports, insulator materials, meals, pastes and any other 
component as j elements. 

In this work we propose a simple model to follow a cooking procedures ac- 
cording with la and H curves in figure [21 This model will help us to design an 
effective furnace for product cooking with a constant energy flow input imple- 
mented by a control system that will have the characteristic of adjusting the 
energy flow input to a constant value Qe = Eq. Here it is important to mention 
that Eq includes all energies we are supply to the furnaces, i.e. Eq involves 
electrical power in case of electrical ovens, or heat coming from combustion in 
case of gas ovens, or solar power in case of solar furnaces or any combination of 
the former for hybrid energy devices. With this assumption we can write 

h-Qs=CT^, (3) 

from Eq. ([ij to describe the temperature of the oven. 

Within this approximation, and according to statements 1 to 5, we will 
estimate the maximum value for Qs by the relation that evaluates heat flow in 
steady state, which can be expressed mathematically as: 

T — T 

Qs = (4) 

where Ta is the environment temperature, and Rt is the total thermal resistance 
from the chamber to the environment. Substitution of Eq. (|4]) in Eq. ([3]) gives 

+ r = Tynaxj (5) 
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where we have defined: 
and 



T = RtCj 



Tmax — EqRt + Ta, 

here r and T^ax are constant. The sohition of this equation is 
Tit) = T„ 





'-t' 




'-t 


^1 — exp 




^ + Tq exp 






T 




T 



(6) 
(7) 

(8) 



The relation ([5]) describes the thermal behavior within the chamber, and by 
selection specific parameters r and Tmax it displays a similar curve to the one 
shown in the figure [2] with the curves lay II. 

To be able to observe the behavior of the difference temperatures between the 
maximum temperatures Tmax and the room temperature Ta, we substitute (O 
in the relation (|S]), and after doing algebraic simplifications we find this relation: 



T{t) ~Ta^ EoRt [l-Exp 
If we use Q = T — Ta and Qmax = 

e{t)^emaxit) [I -Exp 



(9) 



EqRt, then relation © can be written as 
t 



(10) 



As an example we graph the equation ^ varying r and considering Tmax 
constant, and viceversa, we have also considered that the furnace cooks products 
at a temperature of To — 1000°C, and the results are shown in the figures[3]and 
4, which when compared with the normal thermal performance of the curves la 
and // in the figure [T] show good agreement. 

To model the maturation curve /{, is required to supply a constant energy 
Eqj , that compensates the energy dissipated to the environment. For this 
reason to obtain T{t) = To during the time t in the interval tc < t < tc + tm, 
is required: 

To-Ta 



E, 



OJ 



Rt 



(11) 



Finally to model the cooling curve Ic in the interval tc + tm ^ t < tc + tm + te, 
we must consider Ef^ — in the balance equation ([3]), in which, we must also 
substitute the relation (U) , therefore we obtain: 



T-Ta 
Rt 



~dt' 



the relation ([T^ has the solution: 

T^Ta + iTo- Ta) Exp (^^£±^^ Exp (^-^) , 



(12) 



(13) 
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Figure 4: Theoretical performance of the furnace, considering constant Tmax = 
1200°C, and varying r. The hues: 1. t = 0.5/ir; 2. r = Ihr; 3. r = 2hr; 4. 
T = 3hr; A. t = Ahr. 
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or the relation: 



e = QoExp 



~^ 



, -t 

Exp I — 



(14) 



where Qo = To-Ta. 

Here we need to emphasize that E'o.i can be provided to the furnace using 
different energy sources, it can be supplied by electricity, or heat from combus- 
tion, or radiation, or any combination of them and some other energy source. 

Once we have modeled the cooking curves, la, It and Ic', and also la along 
with //, and defining the way of joining the different energy sources, we will 
use the development of furnace sizing. 

Now we calculate the change of energy accumulated in the furnace through 
the relation: 



AEac = CT*0{t) = CT* EoRt [I - Exp 

= tEo ( 1 - Exp 



T 



(15) 



(16) 



(17) 



We can observe that according to these results, the energy dissipate to the room 
is an increasing function on time, and the accumulated energy is bounded, which 
is in perfect compatibility with the physical processes of the furnaces operation. 



The total energy that enters the furnace according to the model is: 

Et ~ Eot, 

therefore the energy dissipate to the room is: 

Qs=Et- Eac = ^0 ( i - T (l - Exp 



2.2 Furnace design. 

We use the previous ideas to design a furnace, according to the figure [2 with a 
rectangular form, we must know, 

• The maximum operating furnace temperature Tq. 

• Cooking time ic- 

• Maturing time tm- 

• The temperature T^. 

With this information, we determine the furnace parameters: thickness of in- 
sulator i5, Eo.i, Eqj, Qmaxi Rt, Ct, t and Tmax to can describe the thermal 
behavior of the oven. 
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Index j 


Material 


p [kg /nv^) 


Cp,, ( J/fcg°C) 








1 


Cooking material 


Pi 


Cp,l 




mi 


Wi = 1 


2 


Heat-resistant 


92 


Cp,2 




7712 


W2 = I 


3 


Thermal insulator 


P-i 


Cp,3 




7773 


W3 = 0.55 



Table 1: Defined parameter for the evaluation of Ct 



The parameter r determines the speed of temperature evolution, and while 
heating or cooking products we can observe two different types of thermal evo- 
lution, as it can be seen in the figure 12.11 a slow temperature curve where r has 
a value in the interval: 

-tc < T < -t,, (18) 
and a fast temperature curve when r has a value on the interval: 

Y^t,<T<-t,. (19) 

Rt and Ct are the parameters that control the thermal behavior, for this, 
it must be adjusted to fit the relation: 

Ct = ^ = f^, (20) 

where a is the selected in the interval < a < | and tc from relations (|18p 
and p9)) . We must also take into account the relations ([2]) and ((20)) to fulfill 
the requirement: 

5Z CpjmjWj = ^ Cp^jPjVjWj = (21) 

3 j 

where the material j has the density pj , mass nij and the volume Vj . We will 
define j materials and their properties according to the table [TJ 

Also we can determine Omax by considering the value of tq = ate and that 
the cooking temperature Qo is reached at the time t = tc, after, when replacing 
in the relation (fTOl) . we get: 



1 - E xp I — 
V a 



(22) 



where Qo — To — Ta- from the relation when doing some algebra we get 
to: ^ 

©max = "ir = / -1 M 1 (23) 
[1-Sxp( — )J 

It is also important to determine T,-aax — ^max + 7a and if we replace the 
value of Qmax from the relation ()23p we obtain: 



[l-^a;p(— )J 
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To determine the energy value for the curve la we know that 



Eq — EqA — 



Ri 



This relation is the adequate energy flow for the furnace operation at a constant 
energy, and we must equal it to the relation of the energy flow of the sources, 

The energy Eqj needed to preserve the temperature during the time interval 
tc < i < ic + is evaluated while considering the relation pil) and the cooking 
conditions, therefore we obtain: 



{To - To) 
Rt 



(25) 



Finally, the temperature T{t) of the furnace, in the interval < t < tc, 
according to the model will have the relation: 



T{t) = Ta 



1 — Exp 



-t 



atr 



[l-Exp{-^)] 
to evaluate 9(t) we can do it through the relation 



TaExp 



-t 



(26) 



[^-Exp{^)] 



1 - Exp 



atr 



(27) 



In the interval t^ < t < t^ + t„i, will have the constant temperature T{t) = 
Tq. In this interval, the cooling temperature will evolve following the relation: 



T(<) = Ta + (To - Ta) Exp 



~f~ 

To 



Exp 



-t 
atr 



or because of the temperature 9 like: 



Q(t) = QoExp 



H~ ^r; 



To 



Exp 



atr 



(28) 



(29) 



Once all the parameters required to characterize a square furnace have been 
met, we will present conclusions. 



3 Oven sizing 

Here we give some examples using different geometries as rectangular, cylindrical 
and semi-spherical. 

We have developed the following design steps to determine all the parameters 
involved in the oven 

1. Let determine t^ and Tq according with the typical cooking curve. 
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Figure 5: Standard oven shapes: a), cylinder., b) semi-spherical, c) rectangular. 
Cooking chamber 1) thermal insulator 2) energy source 3) cooking material 
4)dead mass 5) refractory support. 

2. Let Ti the initial temperature, taken as Ta the ambient temperature. 

3. From T{t) curve describing the expected process we can select r = tq. 

4. The geometrical form of the oven is selected according with practical in- 
formation of the specific process. 

5. Formulas for Rt{5) and Ct{S) are determine according with the geomet- 
rical form and 5 = is a fitting parameter to be determined from 

To = Ct{6)Rt{6). (30) 

6. With ^0 we can obtain Rt{5o) and Ct(<5o)- 

7. With these data we can compare T„iax, Eq, -Eo.m, etc. and the curve T{t) 
with the expected curve. 

8. After some iterations we can get the final geometrical parameters for the 
oven. 

Thus, in order to sizing the oven we need to consider relations 1201 and 1211 to 
determine the maximum load. From these relations we observe that t{6), then 
to determine the adequate thickness we need to propose that 

To = t((5o) = ate. 

where 5o is the selected thickness of the insulator such as with this value the 
corresponding tq to follow the cooking process we select. 

p _ Re,p 

Rt - -J-, 

Zip 

Where Ap is the average area of the wall and i?e,p is the thermal resistance per 
area unit. In the case of a rectangular oven with lengths /i, '3 and thickness 



13 



S we can get: 

" 2 {hh + hh + hh) + 4 (Zi + l2 + h)5+ 126^ ■ ^^^^ 
This expression for the case of a cube is transformed in 

- 6/2 + 12Z^+m2- (32) 
In the case of a cylinder close by to parallel walls, we can write 

1 1 

^"^^"'^ " A^,cu I £ " " 7r[2n,,(r,,e+^)+a^] ^ ' ^^^^ 

^e,p ^e,cil ^e,p ^e,cil 

and in this case the specific Re,cii is 

Ln( ^ ^ ) 2 

2 7r/i;ais 2 7r (ri,cii + d) n. 

where rj^cii is the internal an external radius of the cyhnder and ri^cU is the 
internal radius of the insulator. 

In the case of a semi-spherical oven we have that 

_ 1 _ 1 

-Re,p ^e,esf '^Re,p Re,esf 

where Vi^esf is the radius of the sphere and the thermal resistance can be ob- 
tained as 

"'"^^ An {n^cil + <5) Ti^cil kais 477 {ri,esf +Sf h' 

where Ti esf is the internal radius of the sphere. 

As we need the heat capacities of the insulator for different shapes, in these 
cases we have: 

1. Parallelepiped shape. 

Considering li , I2 and Z3 as the dimensions of the parallepiped and 5 as 
the thickens of the insulators we get 

[85^ + 4(/i +h + h)5'^ + 2{hh + hh + hh)5] 

(37) 

where Cp a^g, pais and 5 are the mean heat capacity, density and thickness of 
the insulator, respectively. Here C^^p is heat capacity of the parallelepide 
oven without considering the insulator term. 
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Component 






mj (kg) 


Cp.jiJ/kg ■ K) 






Product 


Pi 




0.380 


820 


1 


311 


Heat-resistant 


P2 


V2 


2.650 


820 


1 


2.146 


Electrical resistance 


P3 


U3 


0.3 


450 


1 


1.50 


Structure 


P5 


"5 


8.0 


450 


0.05 


180 


Insulator 


P6 


W6 


6.4 


1000 


0.45 


2880 



Table 2: Physical properties of the oven[2] . 

2. Cubic shape. Considering I as the size length, then, 

CT,cub = Cx.cub + Cp^aisPaisgais [^S^ + 12/(5^ + 6PS] (38) 

Here Cx,cub is heat capacity of the cubical oven without considering the 
insulator term. 

3. Cylinder shape. Here L and ri^cU are the length and radius of the cylinder, 
then, 

.aisPaisQais \^ 

(2J3 + (L + An,,u)5'' + 2n,,u{L + 1)5)] 

(39) 

where again Cx,c is heat capacity of the cylindrical oven without consid- 
ering the insulator term. 



4. Semi-sphere. With ri e^/ as radius 



(40) 



where again Cx,e is heat capacity of the semi-spherical oven without con- 
sidering the insulator term. 

With these formulas our description of the time evolution temperature inside 
the oven is complete. 

3.1 Example of a electrical oven 

Using this methodology we have designed an electrical furnace. The correspond- 
ing parameters are given in Table [2j This furnace was use to cook a ceramic 
piece. During the cooking process the temperature was taken and they are plot- 
ted in fig. m as it can be seen, the agreement between model and real data is 
good. 

4 Conclusions 

In this paper we present a simple model to describe the mean temperature 
inside an oven considering global thermal balances. With this model we ana- 
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Figure 6: Comparison between experimental data and themodynamical model 

lyze the thermal behavior under changes of the physical properties of the oven 
components. 

Since different products require different temperature time evohition in the 
cooking process, we present a methodology to design effective ovens, i.e., ovens, 
that from a global point of view, cook following a specific temperature evolution. 
Here is important that the criterion is in the effectiveness instead the efficiency, 
because it is more important following the required temperature than the use 
of the energy in the cooking process. 

We have considered different shapes and present examples of our methodol- 
ogy- 

We test this methodology performing an experiment under specific parame- 
ters and comparison between theoretical and experimental results is good. 
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